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Abstract 

A generalized version of Bertrand's theorem on spherically symmetric curved spaces is 
presented. This result is based on the classification of (3 + l)-dimensional (Lorentzian) 
Bertrand spacetimes, that gives rise to two families of Hamiltonian systems defined 
on certain 3-dimensional (Riemannian) spaces. These two systems are shown to be 
either the Kepler or the oscillator potentials on the corresponding Bertrand spaces, and 
both of them are maximally superintegrable. Afterwards, the relationship between such 
Bertrand Hamiltonians and position-dependent mass systems is explicitly established. 
These results are illustrated through the example of a superintegrable (nonlinear) oscil- 
lator on a Bertrand-Darboux space, whose quantization and physical features are also 
briefly addressed. 



1 Introduction 



Bertrand's theorem, which dates back to the XIX century [I], is a landmark result in 
classical mechanics characterizing the Kepler and harmonic oscillator potentials in terms 
of their qualitative dynamics. More precisely [2], Bertrand's theorem asserts that any 
spherically symmetric natural Hamiltonian system H = \\\>\ 2 + V(|q|) in (a subset of) M 3 
that has a stable circular trajectory passing through each point in its configuration space and 
all whose bounded trajectories are closed is either a harmonic oscillator (V(r) = Ar 2 + B) 
or a Kepler system (V(r) = A/r + B). 

Surprisingly, the classical theorem of Bertrand found a natural extension in and applica- 
tion to general relativity some fifteen years ago thanks to a remarkable paper of Perlick [3]. 
Indeed, the author undertook the classification of all Bertrand spacetimes, which, roughly 
speaking, are spherically symmetric and static spacetimes whose timelike geodesies satisfy 
properties analogous to those of the trajectories of the harmonic oscillator or Kepler sys- 
tems. The connection between timelike geodesies in spacetime and the trajectories of a 
classical Hamiltonian system is that, if one writes the Lorentzian metric as 

r l = g ij (q)dq i dqi -V-Xq) dt 2 , 

where g is a Riemannian metric on a 3-manifold, the timelike geodesies in spacetime are 
naturally related to the trajectories of the Hamiltonian system 

H = g ij (q) PlPj + V(q). 

Perlick's classification of Bertrand spacetimes consisted of two multi-parametric families. 
When one additionally imposes that the spacial part of the metric be Euclidean, Perlick's 
result becomes tantamount to the classical Bertrand's theorem. Moreover, the general case 
includes a number of other systems that have received considerable attention, particularly 
in connection with integrable monopole motion and the existence of generalized Runge- 
Lenz vectors. For instance, Perlick's classification includes spacetimes constructed over the 
3-sphere, the hyperbolic 3-space and the Iwai-Katayama spaces [HE], which generalize the 
Taub-NUT spacetime. 

In this paper we review several recent results related to this problem and we also estab- 
lish a natural connection between curved Bertrand systems and position-dependent mass 
(PDM) Hamiltonians [6]-[20j. We begin by showing, in the next section, that the rather 
complicated families of Bertrand spacetimes admit a strikingly simple physical interpreta- 
tion [21J: they correspond to either an intrinsic oscillator or an intrinsic Kepler system. 
Next we address in section [3] the superintegrability of the associated Hamiltonian system 
and recall the construction of additional integrals of motion of Runge-Lenz type, which 
settled in a satisfactory way a problem with a large body of related literature (cf. [22j [23] 
and references therein). This allows us to state an optimal version of Bertrand's theorem on 
Riemannian manifolds. Furthermore, we present in section 4 the application of the above 
results to PDM Hamiltonians by rewriting the previous Bertrand Hamiltonians in terms 
of a variable mass function. All of these results are explicitly illustrated by discussing in 
section 5 the iV-dimensional (NT)) version of one of the most interesting Bertrand Hamil- 
tonians: the Darboux III Hamiltonian [24]. The associated quantum mechanical problem 
together with some possible physical applications of these new PDM integrable models are 
also considered in the last section. 
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2 Bertrand spacetimes 

To begin with, let us recall Perlick's definition of a Bertrand spacetime [3J. We will consider 
a spherically symmetric, static spacetime (.MxR, 77), where M is a 3-manifold. This ensures 
that the Lorentzian metric r\ can be written as 

dt 2 

T] = h{rf dr 2 + r 2 (d9 2 + sin 2 9 d^ 2 ) 



V{r) ' 

where V is a smooth scalar function and 

g = h(r) 2 dr 2 + r 2 (d9 2 + sin 2 9 dip 2 ) (1) 

defines a Riemannian metric on M. Following Perlick, by a trajectory in spacetime we 
mean the projection of an inextendible timelike geodesic to a (fixed but otherwise arbitrary) 
constant time leaf M x {to}- This terminology is motivated by the fact that a trajectory in 
spacetime actually corresponds to a trajectory (in configuration space) of the Hamiltonian 

in Ai. As customary, g tJ is the inverse matrix to gij, p S T*M is the momentum and 
(j> r ,Pe,Pifi) are the conjugate momenta of the coordinates (r,8,tp). 

Definition 1. The Lorentzian 4-manifold {M x R, rf) is a Bertrand spacetime if: 

1. There is a circular (r = const.) trajectory passing through each point of A4. 

2. The above circular trajectories are stable, that is, any initial condition sufficiently 
close to that of a circular trajectory gives a periodic trajectory. 

Perlick's classification of all Bertrand spacetimes [3] can then be stated as follows: 

Theorem 2 (Perlick). The metric of a Bertrand spacetime can be expressed in exactly one 
of the following forms: 



1. Type I: 

m 2 dr 2 
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n 2 (1 + Kr 2 ) 
2. Type II: 



+ r 2 (d9 2 + sin 2 9dcp 2 ), V = \7r" 2 + K + G. 



2m 2 1 - Dr 2 ± y/(l - Dr 2 ) 2 - Kr 4 

-l 



V = G T r 2 [I - Dr 2 ± - Dr 2 ) 2 - Kr 4 
Here m, n are coprime nonnegative integers and D, G, K are real constants. 
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Notice that in (3j [2T] this result was expressed through the rational number (3 = n/m. 

To provide an interpretation of the above (rather involved) formulas for the Bertrand 
spacetimes, it is convenient to define the concept of harmonic oscillator and Kepler potential 
in any spherically symmetric 3-manifold. To this end, let us use the notation ([1]) for a 
spherically symmetric metric on Ai and consider its associated Laplacian, which we denote 
by A g . It is standard that if u(r) is function on Ai that only depends on the radial 
coordinate, then its Laplacian is also radial and given by 

A ^ d / r 2 du 

9 r 2 h(r) dr \h(r) dr 

Then the symmetric Green function u(r) is obtained as the solution of the equation A g u(r) = 
on .M\{0}, namely 

"to = f dr*. (3) 



As the Kepler potential in 3D Euclidean space is simply the radial Green function u(r) of 
the Laplacian and the harmonic oscillator is its inverse square, it is natural to make the 
following 

Definition 3. The Kepler and the harmonic oscillator potentials in (Ai,g) are respectively 
given by the radial functions 

V K {r) = A l(f a r'- 2 h(r') dr' + B^j , V Q {r) = A 2 ^ jf r'- 2 h{r') dr' + B^j , 

where a,Aj,Bj (j = 1,2) are constants. 

This definition is obviously valid in higher dimensions as well. 

Example 4. Let (Ai, g) be the simply connected 3D space of constant sectional curvature k. 
In this case the metric has the form (fTJ with h[r) 2 = 1/(1 — Kr 2 ), so that the corresponding 
Kepler and harmonic oscillator potentials are 



V K = Vr- 2 - k , V = ( r - 2 - k)- 1 , 

up to additive and multiplicative constants. In terms of the distance function p K to the 
point r = these can be rewritten as 

Vk = Vncot (-v/^Pre) , V = tan 2 (^p K )/K, 

thus reproducing the known prescriptions for the sphere and the hyperbolic space |25[ 126] . 

This readily gives the following interpretation of the Bertrand spacetimes [21] : 



Theorem 5. In a type I (resp. type II) Bertrand spacetime, V is the intrinsic Kepler (resp. 
harmonic oscillator) potential associated with g. 
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3 Bertrand's theorem in 3D curved spaces 



Let us now consider the maximal superintegrability (MS) of the Bertrand Hamiltonians ([2]). 
It is well known that, the Bertrand Hamiltonians being spherically symmetric, to establish 
their MS it suffices to obtain a functionally independent additional integral. As it turns out, 
the most convenient way to obtain this additional integral is as a 'generalized Runge-Lenz 
tensor'. 

Indeed, Bertrand Hamiltonians are somehow similar to the usual harmonic oscillator 
and Kepler systems in Euclidean space. It is classical that the MS of the Kepler system 
can be readily proved using that the Runge-Lenz vector is conserved [27]. In the case of 
the harmonic oscillator, there is no natural way of defining a conserved Runge-Lenz vector, 
but the elements of the symmetric Fradkin matrix Cij = 2u 2 qiqj + piPj are constants of 
the motion and encode the main algebraic properties of the model [28]. A straightforward 
computation shows that the standard Kepler system is obtained from the type I Bertrand 
Hamiltonian by setting K = and n = m = 1, while the harmonic oscillator is the type 
II Hamiltonian with parameters K = D = 0, n = 2 and m = 1. The statement of the 
following key theorem, which proves the MS of Bertrand Hamiltonians [22], is therefore not 
surprising: 

Theorem 6. There exists a (nontrivial) rank-n symmetric tensor field invariant under the 
flow of the Bertrand Hamiltonian, where n is the parameter introduced in theorem \^ 

The proof of this theorem has been fully discussed in [22J. Here we want to stress that 
the parameter n plays a crucial role in the construction of the first integrals. The basic 
observation, which goes back to Fradkin [29], is that any spherically symmetric Hamiltonian 

= f|p| 2 + k r (l f ll) preserves the unit vector field 

cos (p sin ip 2 

a = q H — q x (q x p) , J = p = r ip. 

r rJ 

Of course, this only shows the existence of a local additional integral, which is trivial in 
view of the flow-box theorem, but the point is that this provides a bona fide global first 
integral provided one can express cose/? and sin ip/ J in terms of p and q. In the case of 
the Kepler problem, this readily yields a conserved vector field which is essentially the 
Runge-Lenz vector divided by its norm. In the case of the harmonic oscillator, the above 
quantities are not well defined functions of p and q, but they do define an analogous 
bi-valued conserved vector field. By taking the two-fold symmetrized tensor product of 
this vector field one can remove this indeterminacy, and this procedure yields an invariant 
2-tensor essentially analogous to the conserved matrix Cjj mentioned above. The proof 
of theorem [6] follows this line of thought a bit further. Ultimately, the role n plays in the 
proof merely reflects the properties of the trajectories of the Hamiltonian (studied in [3tl21]). 
which depend crucially on the numbers n and m. It should be noticed that the dependence 
of the constants of motion on the momenta could be extremely complicated, which is the 
reason why only a few among the Bertrand Hamiltonians had previously been identified as 
maximally superintegrable systems. 

To summarize, let us state the complete version of the optimal extension of Bertrand's 
theorem to Riemannian manifolds |22j : 
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Theorem 7. Let H = \g^PiPj + V be a Bertrand Hamiltonian, i.e., a spherically sym- 
metric, natural Hamiltonian system on a Riemannian 3-manifold (M.,g) that has a stable 
circular trajectory passing through each point in its configuration space and whose bounded 
trajectories are all closed. Then the following statements hold: 

1. The metric g and the potential V are of the form given in theorem\^for some coprime 
positive integers n,m. 

2. The potential is the intrinsic Kepler or the harmonic oscillator potential of (A4,g). 

3. H is superintegrable. More precisely, there exists a nontrivial rank n tensor field which 
is invariant and plays the role of the Runge-Lenz vector. 



4 Curved Bertrand systems as classical PDM Hamiltonians 

So far, by starting from Perlick's classification of (3 + 1)D Bertrand spacetimes, we have 
obtained two families of MS Hamiltonians on 3D Riemannian manifolds, which are either 
of Kepler or oscillator type. These results can also be translated into the language of PDM 
systems [6]-|20j. For this purpose, we remark that Bertrand Hamiltonians are initially ex- 
pressed in the form ([2]), that is, H = ^g 13 PiPj + V(r), with the potential V being determined 
by theorem 2. Hence it is necessary to rewrite the Hamiltonian in terms of a variable mass 
M(|q|) in the form 

H =mMi +vm ' (4) 

At this point we stress that the 'radial' Bertrand coordinate r is by no means |q|. Therefore, 
the translation has to be achieved by defining the appropriate change of coordinates r •<-)■ |q| 
in the underlying metric (JTJ) of the Bertrand Hamiltonians, thus giving rise to a conformally 
flat metric: 

g = h{rf dr 2 + r 2 (d0 2 + sin 2 9 dip 2 ) = /(|q|) 2 dq 2 . (5) 
This yields the relations 

r=|q|/(|q|), /(|q|)d|q| = h(r)dr, ^d|q| = ^dr. (6) 

Consequently, the variable mass can thus be related to the conformal factor by setting 

M(|q|)=m /(|q|) 2 = ^, (7) 

where mo is a positive real constant that hereafter we shall fix to 1. In the following we 
present this relationship for each of the two types of Bertrand Hamiltonians. 

4.1 Type I: Bertrand-Kepler Hamiltonians 

In this case, by taking into account that 

h(r) 



n 



Vl + Kr 2 
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and by applying relations ([6]) we obtain that 

m/n 

|q| 



1 + Vl + Kr 2 J ' |q|-("M - K\q\( n / m ) ' 

Therefore the position-dependent mass function is given by 

M(h\) = 5—, (8) 

(|q|" (n/m) -^|q| (n/m) ) q 2 

and the resulting PDM Hamiltonian provided by theorem 2 turns out to be 

1/2 



ff=I (jq|-( n / m ) -K|q|( n / m )) 2 q 2 p 2 + Ai jif + ^ (|qr (n/m) - K|q| (n/m) ) 2 | 



0) 



Example 8. Riemannian spaces of constant sectional curvature k arise in the Bertrand 
Hamiltonians of type I when n = m = 1 and K = —k [21]. Hence the variable mass ([8]) 
and the Hamiltonian Q reduce to 



which is the known Kepler system written in Poincare coordinates [251 130] on the spherical 
(k > 0), hyperbolic (k < 0) and Euclidean (k = 0) spaces. Clearly, we can scale the 
Hamiltonian as 4:H K . The corresponding Runge-Lenz vector in these coordinates can be 
found in EaCSB. 



4.2 Type II: Bertrand-oscillator Hamiltonians 

Now the function h{r) is given by 

mV2 ( l-Dr 2 ± - Dr 2 f - Kr* \ ^ 
h{r) = — { (1 - Br 2 ) 2 - Kr± J ' 

and the relation between the variables |q| and r defined through ([6]) can be obtained as 

|q| = exp < r u(r) — / u(r')dr' > , (11) 



(up to an additive constant coming from the integral) where u(r) is the Green function ([3]). 
For these type II systems u(r) reads as 

u f r ) = (i- Dr 2± J (I - Dr 2 f - kA 1 ' 2 . 

nr V ' / 

Therefore, the explicit general result for M(|q|) = r 2 /q 2 is quite cumbersome and we omit 
it here. However, for some particular cases it adopts a simple form, as it is shown in the 
sequel. 
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Example 9. The three Riemannian spaces of constant curvature k now appear by setting 
n = 2, m = l, K = 0, D = k and by taking the positive sign within the two posibilities '±' 
in (I10p [21]. In this way we find that 



r.2 

V J- rvi 

so that 



h ( r ) = u M ( r ) = 



\f\- HI- 



2q 



1 + Kq 2 ' ' Ml 1 + y/l- — 2 



provided that we have dropped an additive constant In 2 in the integral (|11|) . Then, as ex- 
pected, we obtain the same variable mass as in the previous example and the corresponding 
Bertrand oscillator in Poincare variables is given by 

H K = \(l + K<?f V 2 + A 2 2q2 



(1 -nq 2 ) 2 ' 

The integrals Cu = I{ in the diagonal of the conserved matrix Cy can be found in [25] . 

5 The Darboux III oscillator 

In what follows we focus our attention on a particular system of the family of Bertrand 
Hamiltonians of type II, the so called Darboux III oscillator. The underlying Bertrand space 
is the 3D version of the Darboux surface of type III |31} 132] . for which an iVD spherically 
symmetric generalization was constructed in \33\ 134]. Such a 3D Darboux-Bertrand space 
corresponds to choose the + sign in (I10p and to set n = 2, m = 1, K = D 2 and D = -2A, 
where A is a real parameter. This yields 



, 1 + VI + 4Ar 2 1 + VT+ 4Ar 2 \ 

n(r) = — , uir) = — . 

V ' 2V1 + 4Ar 2 V 2r / 

Therefore the transformations between the radial variables r and |q| turn out to be 

/ / \ 1/2 

/ Vl + 4Ar 2 - 1 \ 



IqlTT+Aq 2 , 



2A 



Then the variable mass function reads M(|q|) = 1 + Aq 2 and the resulting Darboux Hamil- 
tonian is given by 

h p 2 i ^ 2 q 2 no) 

^"2(1 + Aq 2 ) + 2(1 + Aq 2 )' {U) 

where we have written A 2 = u 2 /2. According to section 3, the 3D Hamiltonian H\ is a 
MS system, since it is endowed with a conserved matrix dj (a curved Fradkin tensor). In 
fact, this result can directly be extended to arbitrary dimension N, as it has been proven 
in [SUES]: 



8 



Theorem 10. (i) The Hamiltonian H\ for any dimension N and for any real value 

of X, is endowed with the following constants of motion. 

• (2N — 3) angular momentum integrals: 

C {m) = ^ (qiPj - qjPif , C( m )= ^2 (<HPj ~ QjPi) 2 > ( 13 ) 

l<*<i<m N-m<i<j<N 

where m = 2, . . . , N and C^ N ' = Cnv) • 

• N 2 integrals given by the components of the ND curved Fradkin tensor: 

Cij = PiPj ~ (2AF A (q, P) - u 2 )qiqj, (14) 
where i, j = 1, . . . , N and such that H\ = \ Y^,i=i @u- 

(ii) Each of the three sets {H x , C^}, {H x , C (m) } (m = 2, . . . , N) and {Co} (i = l,...,N) 
is formed by N functionally independent functions in involution. 

(Hi) The set {H\, C^ m \ Cr m \, Co} for m = 2, . . . ,N with a fixed index i is constituted by 
2N — 1 functionally independent functions. 

We remark that the constants of motion (|13p and (|14|) can also be obtained [35J from the 
free Euclidean motion by means of a Stackel transform or coupling constant metamorphosis 
(see [361 EZ] and references therein) . 

It is also worth stressing that although the above statement holds for any real value of 
A, the specific resulting system does depend on such a value, in such a manner that H\ 
comprises, in fact, three different nonlinear physical systems [35J: 

• Nonlinear hyperbolic oscillator. For A > the Darboux space is the complete Rieman- 
nian manifold A4 N = (M. N ,g) with metric gij = (1 + Aq 2 )<5jj. The scalar curvature 
i?(|q|) has a minimum at the origin R(0) = —2XN(N — 1), which coincides with the 
scalar curvature of the A^D hyperbolic space with negative constant sectional curvature 
k = -2A. 

• Nonlinear spherical oscillator. For A < we firstly consider the interior Darboux 
space defined by Ai N = (B\q\ c ,g), where g^ = (1 — |A|q 2 )5jj and -B| q | c denotes the 
ball centered at of radius |q| c = l/-^/|A| (the critical value for the metric and for 
H\). Now R(0) = 2|A|iV(iV — 1) is exactly the scalar curvature of the iVD spherical 
space with positive constant sectional curvature k = +2|A|. 

• Nonlinear exterior potential. For A < we can also consider the exterior Darboux 
space defined by Ai N = (^ N \B^ c , g); this implies to reverse the sign of the metric 
and, therefore, of the Hamiltonian itself, namely, g^ = (|A|q 2 — This system 
can naturally be interpreted as an infinite barrier potential rather that an oscillator 
one. 
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6 Quantization of the Darboux III oscillator 



The quantization problem for Bertrand Hamiltonians arises as a challenging research pro- 
gram, whose first steps in the case of the Darboux III oscillator are summarized as follows. 

Let us consider the quantum Cartesian coordinates and momenta, q, p, with Lie brackets 
and differential representation given by 

d d 2 d 2 
[qi,Pj] =ih6ij, qi = qi, Pi = -ihdi = -ih— , A = 1 h 



dqi ' d 2 qi d 2 q N ' 

Our aim now is to construct the quantum mechanical counterpart of the ND classical 
Hamiltonian (| 1 2 : H\(c[, p) — > H\(c[,p). As it is well known, the crucial point is to obtain 
the quantum analogue of the kinetic term, since there is an order ambiguity in its quanti- 
zation. This task can be faced by applying three different quantization procedures |38j : (i) 
the 'Schrodinger quantization'; (ii) the Laplace-Beltrami quantization (which makes use of 
the Laplace operator on curved spaces); and (hi) a PDM quantization. 

We stress that if we impose that the quantum Hamiltonian H\ keeps the maximal 
superintegrability (that is, the existence of 2N — 2 algebraically independent operators 
that commute with -Ha); then only the Schrodinger quantization yields, in a direct way, 
to fulfill this condition. Nevertheless, the Laplace-Beltrami and PDM quantizations also 
lead to MS quantum Hamiltonians once an additional 'pure' quantum potential term is 
added to the initial quantum Hamiltonian, and such potential terms are related through 
gauge transformations to the Schrodinger quantization. The resulting MS Schrodinger 
quantization of H\ (|12() is characterized as follows [39] (this result is worth to be compared 
with theorem 10). 

Theorem 11. Let H\ be the ND quantum Bertrand- Darboux Hamiltonian given by 

^ = 2aTA«p !+ w^) = WTwj { - R2A + " W) - (15) 

For any real value of A the following statements hold: 

(i) H\ commutes with the following observables: 

• (2N — 3) quantum angular momentum operators, 

Q{m) = ^ ( qi pj - qjpi) 2 , C( m )= ^2 (fcPj ~ QjPi) 2 1 ( 16 ) 

l<i<j<m N—m<i<j<N 

where m = 2 , N and C^ N ' = C'(jv) . 

• A^ 2 operators which form an ND quantum Fradkin tensor, 

Cij = PiPj ~ 2XqiqjH x (q, p) + ^q-Aj, (17) 
where i, j = 1, . . . , N and such that H\ = ^ Qi. 

(ii) Each of the three sets {H x , C^}, {H x , C {m) } (m = 2, . . . , N) and {Co} (i = l,...,N) 
is formed by N algebraically independent commuting observables. 

(Hi) The set {H\, C^ m \ Cr m \, Co} for m = 2, . . . , N with a fixed index i is formed by 2N — 1 



10 



algebraically independent observables. 

(iv) H\ is formally self-adjoint on the Hilbert space L 2 (Ai N ), endowed with the scalar 
product 

(*|*)A= f ^(q)$(q)(l + Aq 2 )dq. 

Jm n 

Clearly, the results of theorem 11 should be adapted to each of the three different 
systems described in the previous section. In particular, we consider here the quantum 
hyperbolic-type Hamiltonian, that has been fully solved in [39j . 

6.1 The nonlinear hyperbolic oscillator: quantum case 

The quantum Hamiltonian H\ with A > has recently been shown to give rise to a new 
exactly solvable quantum model in N dimensions [39] which has both a discrete and a 
continuous spectrum. The discrete spectrum depends on a single principal quantum number 
n = 0, 1, 2 . . . and its eigenvalues are 



2 



= -h 2 \(n + *) \l-Jl+ 5 l. (18) 

Therefore, the degeneracy of this model is exactly the same as in the ND isotropic oscillator 
(whose spectrum is recovered under the limit A — > 0), as it should be expected from the 
beginning due to its MS property. Notice also that the bound states of this system are 
hydrogen-like since 

u? 

lim E n = — , lim (E n+ i - E n ) = 0. 
Therefore H\ has an infinite number of eigenvalues contained in (0, and their only 

2 

accumulation point is |V which is, in turn, the bottom of the continuous spectrum given 

2 

by [|j,oo). The corresponding wave functions can be explicitly found in |39j. 

Finally, a short disgression on possible physical applications of this kind of exactly 
solvable quantum PDM Hamiltonians is in order. Firstly, we recall that the MIC-Kepler 
and oscillator potentials on the 3D sphere have been shown in |10| to be suitable as effective 
models for strong and weak confinement regimes in spherical quantum dots. Secondly, it 
is worth stressing that a parabolic mass function has been proposed in |13t HH] in order to 
describe a one-dimensional quantum well formed by a GaAs/Al^Gai-^As heterostructure. 
In fact, if the concentration x grows in terms of a given spatial coordinate q± as x(qi) = aq\ 
and in this type of material m*(x) = mo(a + bx) [13], then an effective mass function of the 
type m*(qi) = mo (a + baq 2 ) arises. Thus, we have obtained a realistic quantum exactly 
solvable model coming from the Bertrand-oscillator potential on a hyperbolic space with 
non-constant curvature. 
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